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We derive generalised multi-flow hydrodynamic reductions of the nonlocal kinetic 
equation for a soliton gas and investigate their structure. These reductions not only 
provide further insight into the properties of the new kinetic equation but also could 
prove to be representatives of a novel class of integrable systems of hydrodynamic 
type, beyond the conventional semi-Hamiltonian framework. 



1 Introduction 

The generalised soliton-gas kinetic equation represents an integro-differential system [2] 
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Abstract 
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Here /(r/) = f{ri,x,t) is the distribution function and s{ri) = s{rj,x,t) is the associated 
transport velocity. The (given) functions S{r]) and G{ri,fi) do not depend on x and t. 
The function G{r],fi) is assumed to be symmetric, i.e. G{r],fi) = G{fL,ri). 
System (1), (2) with 



S{v) = 4r/ , G{r],fi) = log 



(3) 



was derived in [1] as an infinite-genus thermodynamic limit of the Whitham modulation 
equations associated with the KdV equation, ipt — Qipipx + fxxx = 0, and was shown to 
describe macroscopic dynamics of a soliton gas, a disordered infinite-soliton ensemble of 
finite density [1]. In the KdV context, 77 > is a real- valued spectral parameter and 
the function f{ri,x,t) is the distribution function of solitons over the spectrum so that 
K = f{r])dr] = 0{1) is the spatial density of solitons. If k ^ 1, the first order 
approximation of (2), (3) yields Zakharov's kinetic equation for a dilute gas of KdV 
solitons [11]. The quantity S{ri) in (2) has a natural meaning of the velocity of an 
isolated (free) soliton with the spectral parameter r] and the function G{r],^)/ri is the 
expression for a phase shift of this soliton occurring after its collision with another soliton 
having the spectral parameter fi < rj. Then s{'r],x,t) acquires the meaning of the self- 
consistently defined mean local velocity of solitons with the spectral parameter close to 
7]. A straightforward physical derivation of the kinetic equation (1), (2) for integrable 
systems, based on the original Zakharov [11] phase-shift reasoning was proposed in [3]. 

In recent paper [2], the multi-flow hydrodynamic reductions of the kinetic equation 
(1), (2) were studied using the so-called 'cold-gas' ansatz 



N 



f{rj,x,t) = Y,rMSiv-Vn, (4) 



m=l 



where the 'spectral' components rj^ > r]^~^ > ■ ■ ■ > rj^ > are arbitrary numbers. These 
'isospectral' cold-gas reductions were shown to have the form of systems of hydrodynamic 
conservation laws 

nj = (nV),, z = l,...,iV, (5) 

where the conserved 'densities' = ?7*/(?7*, x, t) and the associated velocities = 
— s(?7*,a;, t) are related algebraically: 

v' = ^i + J2 e'"M'"(t;'" - v') , e"" = e^' . (6) 



Here 



The isospectral cold-gas reductions (5), (6) were proven in [2] to represent integrable 
(semi-Hamiltonian [!)]) linearly degenerate hydrodynamic type systems (see [5], [7]) for 
arbitrary A^, which is a strong indication that the full kinetic equation (1), (2) could 
constitute an integrable system in the sense yet to be explored. 
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The present paper is devoted to a more general multi-flow hydrodynamic approxima- 
tion of the kinetic equation (1), (2), which we derive by considering an ansatz (see, for 
instance, [10]) 

N 

f{rj, x,t) = J2 r t)^iv - V'^i^, t)) (8) 

m=l 

with the 'spectral components' rj'^ = 'r]^{x,t) being (unknown) functions of x and t rather 
than arbitrary constants as in (4). We show that the corresponding iV-flow non-isospectral 
hydrodynamic reductions have the form of 2A^-component hydrodynamic type systems 

ni = (nV),, 4=v%, ^ = l,2,...,iV, (9) 

where the functions u^{x,t), v''{x,t) and r]^{x,t) are related algebraically by the same 
equations (6), (7) provided certain restrictions on the behaviour of the kernel function 
G{ri, fi) for ?7 — )■ /i are satisfied. 

System (9), (5), (6) is not integrable by the standard Tsarev generalized hodograph 
method, because it possesses just Riemann invariants and has double characteristic 
velocities. However, having in mind that this system is obtained as an exact reduction of 
an integrable system (at least for S{r]), G{r],^) defined by (3) — the KdV case), one can 
expect that the multi-flow reductions (9) will be integrable by some modification of the 
generalised hodograph method [9]. This could lead to an extension of the conventional 
notion of an integrable system of hydrodynamic type. We are going to investigate this 
problem in detail in future publications. 



2 Generalised hydrodynamic reductions 
2.1 Evolution equations 

Substituting (8) into (1) we obtain (hereafter we shall be using a shorthand notation 77* 
for ?7*(x, t)) 



d_ 
dt 



fix, mv - v'?j + ^ (^s{r^. X, t) fix, mv - v'?j = 0, 



Differentiating and collecting the terms for 6{r] — rj^) and 5'(?7 — 77*) we obtain 

N N 

^[/f + isiv,x,t)nMv - vl - Et/"^" + siv,x,t)rv:]s'{v - vi = o. (lo) 



n=l n=l 



Here /* = f^{x,t). Evaluating asymptotic behavior of this expression near each point 77* 
we arrive 2N component hydrodynamic type system (cf. (5)) 

fi + {s{r]\x,t)r), = 0, r]l + s{r]\x,t)r]l = 0, n = l,...,N. (11) 

It is instructive to derive the hydrodynamic reduction (11) by a direct calculation. This 
is done by integrating (10) with respect to rj over a small vicinity of each point t] = rj^ 
with the weights 1 and (77 — r/*), respectively. 
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Let us fix x = Xq. Then, assuming the ordering rj^ > t]^~^ > ■ ■ ■ > 77^ > to hold for 
all t in a small vicinity of Xq we introduce closed intervals (Tj = [r/* —Ei, 77* + choosing 
Si > in such a way that in the vicinity of Xq one has rj^{x, t) e cxj iff j = i. 

We now integrate (10) over af. 



N 



N 



n=l 



n=l 



dr] = 0, 



which reduces, after integrating the term with 5'*) by parts, to 



N 



J2[f^ + s{v,x,t)l 



n=l 



ox or] 



dr] = 0. (12) 



Now, integration over cTj immediately leads to the hydrodynamic conservation law: 

fi + {s{r]\x,t)n, = 0, (13) 

which is valid in the small vicinity of xq. If we assume that the above restrictions on the 
behaviour of functions ri^{x,t) hold for any x = Xq G M, t > 0, equation (13) will be valid 
on the entire real line. Setting i = 1,. . . ,N in (12) we immediately obtain the first 
equations in (11). 

To derive the second set of equations (11) we multiply (10) by {rj — 77*) and integrate 
over aj to get 



TV 



j2[f- + {s{v,x,t)nMv-vn{v 



n=l 



N 



n=l 



dr] 



dr] = 0. 



(14) 



If j = i, the first integral in (14) vanishes, while the second one, after integrating by parts 
and utilising the fact that each interval Uj contains only its "own" value 77*, yields 



d_ 
dr] 



{{V - V") [fit + sir]. X, t)fr]l]) 6{r]- r]')dr] 



[fX + s{r],x^)f%] + {r]- ry^)M|M/^^^ ) S{r] - r]^)dr] = 0. 



(15) 



Evaluating the integral in (15) we get 

r]i + s{r]\x,t)r]i = 0, i = 1, . . . , N . 



(16) 



It is not difficult to see that, if j 7^ i, we recover equations (13). Thus, the compatibility 
of the non-isospectral ansatz (8) with the kinetic equation (1), (2) imposes restrictions 
(16) on the functions r]\x,t). 
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Generally, 2 A^- component hydrodynamic type system (11) possesses conservation 
laws 

dtiip,{r]')r) + {s{r]\x,t)ip,i7]')f% = 0, 

where ^ii'rf) are arbitrary functions of a single variable. It is convenient to choose </?i(?7*) = 
rf so that (11) reduces to (cf. (5)) 

Mi = (nV),, r/; = t;Y, 2 = l,...,iV, (17) 

where = vt f \ v"^ = —s{r]\x,t). 



2.2 Closure relations 

The closure relations connecting the field variables u\ v\ and r/* in (17) are obtained by 
substituting the same ansatz (8) into the integral equation (2). Since we are going to use 
the variables instead of /*, we slightly modify ansatz (8) as follows 

N 



i=l 



Substitution of (18) into (2) yields 



TV 

s{rj, X, t) = S{rj) + ^'"^^^^|^[^(^"^' ^> t) - 'iv, ^, t)] , (19) 



U 

m=l 

As in [2], we introduce (see (7)) 



= .^k. (20) 



There is an important point to be made. In the linearly degenerate reductions (5), (6) as- 
sociated with the isospectral ansatz (4) involving arbitrary constants rj^, the dependencies 
of ^* = — S'(?7*) and e^^ on the relevant components of the vector rj = {77^, 77^, ... , rj^} are 
not important from the viewpoint of integrability — these only provide the connection 
with the original nonlocal equation (2) — see [2]. However, under the generalised ansatz 
(8), ?7*s become dependent variables, 77* = ri'^{x,t), so the aforementioned dependencies 
become essential for the structure of the corresponding hydrodynamic reductions. 

Now we pass to the limit as 77 — )■ 77* in (19) to obtain, assuming lim^^^^ s{r], x, t) = — f * 
(continuity), 



= ^e™«-(^^ - ^-) - ^(7^^) + ^ hm G(7/, v^)isiv, x, t) - siv\ x, t)). (21) 



If the limit 



lim G(7/, r]'){s{ri, x, t) - s{r]\ x, t)) (22) 

r/— >-?7' 
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exists then (21) becomes 

= ^e™«"^(^^ - vn - Siv') + 9^in, V, ri), (23) 

where 

gi{u,v,r]) = -—\imG{ri,7]')[s{ri,x,t) - s{r]\x,t)]. 

The existence of the hmit (22) imphes that the function Git], fi) has at most a simple pole 
singularity on the diagonal fJ^ = rj. 

If the limit (22) vanishes for alH = 1, . . . , (which happens if G{ri, /i) either vanishes 
itself or has a singularity weaker than a simple pole as /x — )■ r/) then Qi = and equation 
(23) reduces to the closure conditions (6), (7) obtained for the isospectral cold-gas reduc- 
tion. Below we restrict our consideration just to this, most important, case, which arises, 
in particular, in the case of the kinetic equation for the KdV solitons [1], when the kernel 
function G{ri,fi) has only logarithmic singularity on the diagonal - see (3). 

In conclusion of this section we note that nonexistence of the limit (22) for some 
given G(/i,?7) signifies incompatibility of the delta function ansatz (18) with the integral 
equation (2) for that particular kernel G{fi,ri). 



3 The structure of generalised multi-flow hydrody- 
namic reductions 

Motivated by the results of [2] for the isospectral cold-gas hydrodynamic reductions (5), 
(6) we introduce a symmetric matrix e = [e'^^]NxN with the off-diagonal elements e^^{rf) 
defined by (20) and the diagonal elements e^^ being some new field variables r'^(u, ry). 
Theorem 1. ([2]): Algebraic system (6) admits the parametric solution 



N ^ N 

m' = 5^/3™, v'=-Y,UPm^, (24) 



m=l m=l 



where symmetric functions (3ik{r,r]) are the elements of the matrix (3 = [(3mn]NxN such 
that (3e = —1. 

Proof: We replace (6) by an equivalent system 



TV 



6 + X]e™M"*(w™-t;0. (25) 



m=l 



(note that summation in (25) goes over all m including m = i (cf. (6)). Then (25) can 
be re-written in the form 



TV \ TV 



m=l / m=l 
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Substituting (24) into the above formula we obtain 

(N N \ N N 

1 + = + $^$^en^/3^fce'=\ (26) 

m=lfc=l / m=lk=l 

Taking into account /3e = —1, one can see that expressions at both sides of (26) vanish 
independently. Thus (26) is an identity, hence the parametric representation (24) is 
consistent with system (6). The Theorem is proved. 

Corollary: The field variables r'^(u, rj) are rational functions of the conserved densi- 
ties M™, namely, 



r^ = ^ 1 + ^u'^e'"'=(r7) I , A; = 1, 2, iV. (27) 

Indeed, multiplying both sides of the first relationship in (24) by e^^ and performing 
summation over i we obtain: 

N N N 



m=l Tn=ln=l 

Thus, 

which immediately yields (27). 

Theorem 2: Under parametrization (24) the 2N -component hydrodynamic type sys- 
tem (17), (6), (7) reduces to a quasi- diaginal form: 

4 = v'vl ^ = l,...,iV; (28) 

a,nk 

_ ..k^k , ^ I Y-.,n(^,n_^,fc)___^. I^fc^ k = l,...,N. (29) 
%^k 



i + -k 



Proof: The evolution equations (28) for 77* (x, t) are the same as in (17) so we need only 
to derive equations (29) for r''{x, t), k = 1, . . . ,N. Substituting parametric representation 
(24) into the conservation laws (17) we obtain 



\m=l J \m=l J 



Multiplying both sides by e*'^, performing summation over the repeated index i and using 
the relationship /3e = —1, one arrives at the equation 



N N N N 



Y^T^f^rrndte^' = $^$^em/3m^5xe*' - d,^k. (30) 



A simple but not entirely trivial calculation using the following obvious property of the 
matrix e(r, 77): 

~T{ ^nk^ks i 

and the evolution equations (28) for rj^, shows that (30) reduces to 



.k^k 



s=l \n=lm=l / 



k = i,...,N. (31) 



-'mi U V , 



Then taking into account (see (24)) 

N N 
m=l 771=1 

equations (31) reduce to the form (29). The Theorem is proved. 

Eliminating from (24), we obtain expressions relating v to 77 and r: 

N 

mi 

v\r,r,) = "-^ , (32) 

m=l 

while n*(r,77) are given by the first equation in (24). Now, system (28), (29) is closed. 

For the isospectral case, when rj^ , i = 1,2, N are constants so dft]^ = d^Vj^ = and 
we recover from (29) the Riemann invariant representation 

rl = v\r)ri (33) 

of system (5), (6) obtained in [2] with the use of the linear degeneracy of the isospectral 
cold-gas hydrodynamic reductions. As a matter of fact, one can see now that the Riemann 
invariant equations (33) could be readily obtained directly from system (5), (6) by the 
substitution into it of the parametric solution (24). 

Thus, the 2A^-component hydrodynamic reduction (17) admits parametrization (24) 
resolving algebraic system (6) and reducing the evolution equations to the form (28), 
(29). System (28), (29) has double characteristic velocities f'^( r, 77) (32). However, in 
the general case, just N functions 7]^{x,t) are Riemann invariants (i.e. only a half of 
the complete hydrodynamic system (17) can be written in diagonal form), while the field 
variables r^{x,t) become Riemann invariants only if the corresponding r]^ =const. 

In conclusion of this Section we note that linear degeneracy of system (33) proved in [2] 
implies that (9i;*(r)/(9r* = for all i = 1, . . . A^. The latter property clearly remains valid 
for the characteristic velocities u*(r,r/) of the generalised hydrodynamic reductions (28), 
(29), however, now this is no longer associated with the notion of linear degeneracy of a 
hydrodynamic type system in the classical sense [6], [8] since are no longer Riemann 
invariants and also dv^{v,rji) / drj^ 7^ 0. 
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4 Conclusion 



In this paper, we have derived the generahsed hydrodynamic reductions of the nonlocal 
kinetic equation for a soliton gas (1), (2) by considering the non-isospectral multi-flow 
ansatz (8) for the distribution function. These new reductions have turned out to have 
rather unusual structure which we have revealed by using the parametric solution (24) 
to the algebraic closure conditions (6), (7). More precisely, the non-isospectral A^-fiow 
hydrodynamic reductions of the kinetic equation are shown to represent 2A^-component 
half-diagonal systems of hydrodynamic type (28), (29) with Riemann invariants and 
double characteristic velocities. The feature that makes the derived reductions deserving 
special attention is that, while they are clearly not integrable by Tsarev's generalised 
hodograph transform method [9], they could still prove to be integrable in some new sense 
yet to be understood. Indeed, having in mind that system (28), (29) can be derived as a 
generalised hydrodynamic reduction of the kinetic equation associated with an integrable 
equation (e.g. with the KdV equation — for S{ri), G{rj, n) defined by (3)), one can expect 
that this reduction will be integrable by some nontrivial extension of the generalised 
hodograph method. 

As a by-product of our calculations we recover the Riemann invariant structure of 
the isospectral cold-gas reductions (5), (6), (7) studied in [2]. We note that our present 
compact derivation, unlike that in [2], does not make any use of the linear degeneracy 
property of the reductions under study. 
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